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Abstract 



Extending our previous work in the strictly parabolic case, we 
show that a linearly unstable Lax-type viscous shock solution of a 
general quasilinear hyperbolic-parabolic system of conservation laws 
possesses a translation-invariant center stable manifold within which 
it is nonlinearly orbitally stable with respect to small H per- 
turbations, converging time- asymptotically to a translate of the un- 
perturbed wave. That is, for a shock with p unstable eigenvalues, we 
establish conditional stability on a codimension-p manifold of initial 
data, with sharp rates of decay in all L^. For p = 0, we recover the re- 
sult of unconditional stability obtained by Mascia and Zumbrun. The 
main new difficulty in the hyperbolic-parabolic case is to construct an 
invariant manifold in the absence of parabolic smoothing. 
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1 Introduction 

In this paper, extending our previous work in the semilinear and quasilinear 
parabolic case [Z5, Z6], we study for general quasilinear hyperbolic-parabolic 
systems of conservation laws including the equations of compressible gas dy- 
namics and MHD the conditional stability and existence of a center stable 
manifold of a hnearly unstable viscous Lax shock. This is part of a larger 
program initiated in [TZl, TZ2, TZ3, TZ4, SS, BeSZ, Z5] going beyond 
simple stability analysis to study nontrivial dynamics and bifurcation, and 
associated physical phenomena, of perturbed viscous shock waves in the pres- 
ence of linear instability. As discussed for example in [AMPZ, GZ, Z7], such 
conditionally stable shock waves can play an important role in asymptotic be- 
havior as metastable states, and their center stable manifolds as separatrices 
bounding the basins of attraction of nearby stable solutions. 

The main issue in the present case is to construct a (translation-invariant) 
center stable manifold about a standing viscous shock wave; once this is done, 
the conditional stability analysis follows in straightforward fashion by a com- 
bination of the semilinear argument of [Z5] and the unconditional stability 
analysis of [MaZ2, MaZ3, Z2] in the hyperbolic-parabolic case. As discussed 
in [Z6], the technical difficulty in constructing the center stable manifold 
for quasilinear equations is an apparent loss of regularity in the usual fixed 
point argument by which the center stable manifold is constructed. See also 
the related discussion of [LPSl, LPS2]. We accomplish this by the method 
introduced in [Z6] in the quasilinear parabolic case, combining an implicit 
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fixed-point scheme with suitable time-weighted if^-energy estimates. For 

related energy estimates, see [TZ3]. 
Consider a viscous shock solution 

(1.1) U(x, t) = U(x - St), lim U{z) = U±, 

2:— >±oo 

of a hyperbolic-parabolic system of conservation laws 

(1.2) Ut + FiU),^iBiU)U,%, 

X eR, U, F eW, B e M"''". Profile U satisfies the travehng-wave ODE 

(1.3) B{U)U' ^ F{U)- F{U^)- s{U -U_). 
Denote A{U) = Fu{U), 

(1.4) B± := lim B{z) = B{U±), A± := lim A{z) = Fu{U±). 

z-*±oo z— >±oo 

Following [TZ3, Z2, Z3], we make the structural assumptions: 

(Al) = ( ^ - ( t )' ^ ^ ( ° )' ^ "°"''"Sular, 

where U e Ui E R"-^ U2 G M^ and b G M^'^^'; moreover, the Ui- 
coordinate Fi{U) of F is linear in U (strong block structure). 

(A2) There exists a smooth, positive definite matrix A^{U), without loss 
of generality block-diagonal, such that A^^Au is symmetric, ^22^ positive 
definite but not necessarily symmetric, and (A°A)-|- is symmetric. 

(A3) No eigenvector of A± lies in KerS± (genuine coupling [Kaw, KSh]). 

To (Al)-(A3), we add the following more detailed hypotheses. Here and 
elsewhere, (t{M) denotes the spectrum of a matrix or linear operator M. 

(HO) F,B eC\ k>A. 

(HI) a{Aii) real, constant multiplicity, with a{Aii) < s or a{Aii) > s. 
(H2) o'{A±) real, simple, and different from s. 

Conditions (A1)-(A3) are a slightly strengthened version of the corre- 
sponding hypotheses of [MaZ4, Z2, Z3] for general systems with "real", or 
partially parabolic viscosity, the difference lying in the strengthened block 
structure condition (Al): in particular, the assumed linearity of the Ui equa- 
tion. Conditions (H0)-(H2) are the same as those in [MaZ4, Z2, Z3]. The 
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class of equations satisfying our assumptions, though not complete, is suffi- 
ciently broad to include many models of physical interest, in particular com- 
pressible Navier-Stokes equations and the equations of compressible magne- 
tohydro dynamics (MHD), expressed in Lagrangian coordinates, with either 
ideal or "real" van der Waals-type equation of state as described in Appendix 
A. See [TZ3] for further discussion/generahzation. 

Remark 1.1. Conditions (H1)-(H2) imply that U± are nonhyperbolic rest 
points of ODE (1.3) expressed in terms of the C/2-coordinate, whence, by 
standard ODE theory, 

(1.5) |a;(^- t/±)(a;)| < Ce-'^l^l, < r < /c + 1, 
for X ^ 0, some 77, C > 0; in particular, |C/'(x)| < Ce"'''^'. 

By the change of coordinates x ^ x — st, we may assume without loss 
of generality s = as we shall do from now on. Linearizing (1.2) about the 
(now stationary) solution U = U yields linearized equations 

(1.6) Ut = LU := -{AU), + (BU,),, 

(1.7) B{x) := B{U{x)), A{x)V := dF{U{x))V - {dB{U{x))V)U^, 

for which the generator L possesses [He. S. ZH] both a translational zero- 
eigenvalue and essential spectrum tangent at zero to the imaginary axis. 

Our first main result is the existence of a translation-invariant center 
stable manifold about U. Define the mixed norm 

(1.8) |C/|^i,2 := ^J\U,\l, + \U2\l, 
suggested by the hyperbolic-parabolic structure (Al). 

Theorem 1.2. Under assumptions (Al)-(A3), (H0)-(H2), there exists in an 
H^'^ neighborhood of the set of translates of u a codimension-p translation 
invariant Lipschitz (with respect to H^'"^) center stable manifold /Acs, tangent 
to quadratic order at U to the center stable subspace I^cs of L in the sense 
that 

(1.9) |n„(t/ - u)\h^,. < c\Ucs{u - u)\l^,. 



4 



for U e Aics where Ties and n„ denote the center- stable and stable eigenpro- 
jections of L, that is (locally) invariant under the forward time- evolution of 
(1.2) and contains all solutions that remain bounded and sufficiently close to 
a translate of tJ in forward time, where p is the (necessarily finite) number 
of unstable, i.e., positive real part, eigenvalues of L. 

Next, specializing a bit further, we add to (H0)-(H3) the additional hy- 
pothesis that M be a Lax-type shock: 

(H3) The dimensions of the unstable subspace of A_ and the stable 
subspace of sum to n + 1. 

We assume further the following spectral genericity conditions. 

(Dl) L has no nonzero imaginary eigenvalues. 

(D2) The orbit U{-) is a transversal connection of the associated standing 
wave equation (1.3). 

(D3) The associated inviscid shock (C/_, U+) is hyperbolically stable, i.e., 

(1.10) det(rr, . . . , rp_i, r++i, . . . , r+, {U+ -U.))^ 0, 

where r]~, . . . rp_^ denote eigenvectors of associated with negative eigen- 
values and rp_j_^,...r+ denote eigenvectors of associated with positive 

cigcnvahics. 

As discussed in [ZH, MaZl], (D2)-(D3) correspond in the absence of a spec- 
tral gap to a generalized notion of simplicity of the embedded eigenvalue 
A = of L. Thus, (D1)-(D3) together correspond to the assumption that 
there are no additional (usual or generalized) eigenvalues on the imaginary 
axis other than the transational eigenvalue at A = 0; that is, the shock is not 
in transition between different degrees of stability, but has stability properties 
that are insensitive to small variations in parameters. 

With these assumptions, we obtain our second main result characterizing 
the stability properties of u. In the case p = 0, this reduces to the nonlinear 
orbital stabihty result estabhshed in [MaZl, MaZ2, MaZ3, Z2]. 

Theorem 1.3. Under (A1)-(A3), (H0)-(H3), and (D1)-(D3), u is nonlin- 
early orbitally stable under sufficiently small perturbations in L^nH^ lying on 

the codimension p center stable manifold M.cs of u and its translates, where 
p is the number of unstable eigenvalues of L, in the sense that, for some ol{-), 
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all LP, 



(1.11) 



\U{x,t) - U{x 
\U{x,t)-U{x 



«(t))Up < C(l + t)-^(i-F)|[/(a;,0) - [/(a;)Uin//3, 

a{t))\HB<Cil + tr-^\Uix,0)-Uix)\Linm, 

ait) < C(l + t)-5|f/(x, 0) - U{x)\Linm, 
a{t)<C\U{x,Q)-U{x)\L^nH^. 



Moreover, it is orbitally unstable with respect to small H^'^ perturbations not 

lying in A4, in the sense that the corresponding solution leaves a fixed-radius 
neighborhood of the set of translates of U in finite time. 

Remark 1.4. It is straightforward, combining the pointwise argument of 
[RZ] with the observations of [Z6] in the strictly paraboUc case, to extend 
Theorem 1.3 to the case of nonclassical under- or overcompressive shocks 
with (D1)-(D3) suitably modified as described in [Z6]. This gives at the 
same time new information even in the Lax case, including convergence of 
the phase a to a limit q;(-|-oo) at rate (1 + 1)'^^"^, with \di\ < C(l + 1)^^. 

1.1 Discussion and open problems 

For results of a similar nature in the context of the nonlinear Schrodinger 
equation, but obtained by rather different techniques, we refer the reader 
to [Sc]. The method of [Sc] consists of constructing the stable manifold to 
a known center manifold using detailed linearized and nonlinear estimates 
similar to those used to prove stabihty in the spectrally stable case. That is, 
the existence and stability arguments are carried out at the same time. 

This makes possible a somewhat finer analysis than the one we carry 
out here (in keeping with the somewhat weaker stability properties for the 
nonlinear Schrodinger equation, for which the linearized operator about the 
wave possesses essential spectrum on the whole of the imaginary axis lying 
outside a finite ball), but also a somewhat less general one; in particular, it is 
required that the linearized operator possess no pure imaginary eigenvalues 
other than zero. Note that we make no assumptions on the pure imaginary 
spectrum of the linearized operator in our construction of the center stable 
manifold, but only in the proof of conditional stability. 

An interesting issue pointed out in [TZ3] is that the strong block structure 
assumptions (A1)-(A3) used to obtain the key variational energy estimates 
on which our arguments are based hold for the equations of gas dynamics 
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or MHD written in Lagrangian coordinates, but not in Eulerian coordinates, 
and indeed the energy estimates themselves do not hold in that case; see Ap- 
pendix A, [TZ3] for further discussion. On the other hand, the Lagrangian 
and Eulerian formulations are equivalent by a uniquely specified (up to trans- 
lation, which may be factored out) transformation, and so the results do hold 
in Eulerian coordinates. It is an interesting open question whether the re- 
sults of this paper hold for general symmetric hyperbolic-parabolic systems 
as defined in [MaZ4, Z2, Z3], which include among other examples both La- 
grangian and Eulerian formulations. 



2 Existence of Center Stable Manifold 

Defining the perturbation variable V :— U — U, we obtain after a brief 
computation the nonlinear perturbation equations 

(2.1) Vt-LV = N{V), 

where L as in (1.6) denotes the linearized operator about the wave and 

N = N'{V)x is a quadratic order residual, with (recalling (Al)) 

(2.2) 

Af{V) := (^B{lf + V){U + V), - B{U)U, - B{U)V, - {dB{U)V)U,^ 
- (^F(U + V)- F{U) - dF{U)V^ 

We seek to construct a Lipschitz i^^'^-local center stable manifold about the 
equilibrium = 0, that is, a locally invariant Lipschitz manifold tangent at 
V = to the center stable subspace Scs, that is quadratic-order tangent in 
the sense that |nu\/|^i,2 < C|ncsl^|^i,2. 



2.1 Preliminary estimates 

Conditions (A2)-(A3) imply that {i^A - i'^B)^ < for some 9^ > 

0, for all ^ G M [KSh], which in turn implies by standard spectral theory [He] 
that the essential spectrum of L lies entirely in the neutrally stable complex 
half-plane 3?A < 0, and the unstable half-plane > either contains finitely 
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many eigenvalues of finite multiplicity, or else consists entirely of eigenvalues, 
as is easily shown (by energy estimates, for example) to be impossible. Thus, 
there is a well-defined unstable subspace S„ consisting of the direct sum of 
the finitely many generalized eigenfunctions in the positive half-plane, and an 
associated unstable eigenprojection U^. Likewise, there is a complementary 
center-stable subspace determined by the range of the complementary 
projection U^s := Id — n„. 

Proposition 2.1 ([MaZ3, Z2, Z3]). Under assumptions (Al)-(A3), (HO)- 
(H2), L generates a C° semigroup e^* satisfying 



(2.3) , . 



for some (3 > 0, and for all uj > 0, for all t > 0. 

Proof. Straightforward energy estimates yield resolvent bounds verifying the 
C°-semigroup property; see [MaZ3, Z2] and especially [Z3], Prop. 3.6. The 
same estimates yield 

\{X-L)-%.^L^<C{v) 

for A > 77 and |A| sufficiently large, r] > arbitrary, whereupon (2.3) (i) follows 
by Priiss' theorem restricted to the center-stable subspace [Pr, ProK, KS]; 
see [Z3], Appendix A for further discussion, in particular Rmk. 6.23. Bound 
(2.3) (ii) follows by equivalence of norms for finite-dimensional spaces and 
standard finite-dimensional ODE estimates. □ 

Introducing a C°° cutoff function 

p{x) = 




and recalling mixed-norm definition (1.8), let 

1^1 



N\V) ■.^p(yf^)N{V). 



Lemma 2.2. Assuming (A1)-(A3), (H0)-(H2), the map : H^^^ is 
C'''^^ and its Lipschitz norm with respect to V is 0{5) as 5 — > 0. Moreover, 

(2.4) \N\V)\l2<C\V\1^,,. 
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Proof. The norm in H^'^ is a quadratic form, hence the map 

is smooth, and A^^ is as regular as N. By Moser's inequahty, 
\N{V)\l2 < C(\V\hAV\loo + \d,V2\m\VM + IdML^lVlml) < C\V\l,,, 

and, similarly, \N{V^) - N{V^))\l2 < C\V^ - V^\m,2\supj^^^2C\V^m,2\ and 
thus |dA^|i/i,2^i2 < C|y |i/i,2 so long as |y |i7i,2 remains bounded, in particular 
for \V\ma < 5. Thus, \N\V)\l2 < \N{V)\l2 < C\V\]j,,, for \V\m-2 < 5, 
while N\V) = for \V\m,2 > S, verifying (2.4). The Lipschitz bound 
follows, likewise, by 



N\V^) - N\V'')\l2 < 



P 



V^\m,2 



mv')\,2 



+ 



\V'\m.2 



< 3\V^ - 



Hi.2 



X sup ' + 



sup |(iAr(V^)|^i.2^i2). 



□ 



Corollary 2.3. Under assumptions (Al)-(A3), (H0)-(H2), 



(2.5) 



for some (3 > 0, and for all cu > 0, for all t >0, with Lipschitz bounds 



(2.6) 



\e-^^U^dN^\m.2^m,2 < CJe'^^ 



2.2 Fixed-point iteration scheme 

Applying projections Hj, j — cs, u to the truncated equation 

(2.7) Vt- LV ^ N\V), 
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we obtain using the variation of constants formula equations 

UjV{t) = e^^^-'^^'^^Vitoj) + I e^^'-''^IijN\V{s)) ds, 

Jtoj 

j — cs,u, so long as the solution V exists, with toj arbitrary. Assuming 
growth of at most |y(i)|ij2 < Ce^* in positive time, we find for j — u 
using bounds (2.3) (ii) and (2.6) (ii) that, as to,« +00, the first term 

e^*^*~*'''"^n„\/(to,M) converges to zero while the second, integral term converges 
to J^^°° e^^^-'^UuN\V{s)) ds, so that, denoting w := UcsV, z := n„1/, we 
have 

/+00 
e^^^-'^U^N^{{w + z){s)) ds. 

Likewise, choosing to,cs = 0, we have 

(2.9) w(t) = e'^'UM + f e^(*-^)n,,7V^((«; + z){s)) ds, 

Jo 

wq := IlcsV{0). On the other hand, we find from the original differential 
equation projected onto the center stable component that w satisfies the 
Cauchy problem 

(2.10) wt - Lw = UcsN\w + z) 

with initial data wq = IlcsV{0) given at t = 0. We shall use these two repre- 
sentations together to obtain optimal estimates, the first for decay, through 
standard linear semigroup estimates, and the second for regularity, through 
the nonlinear damping estimates (2.14)-(2.15) and (2.18)-(2.19) below. 

Viewing (2.9), or alternatively (2.10), as determining w = yV{z,Wo) as a 
function oi z, we seek 2; as a solution of the fixed-point equation 

(2.11) z = f{z. Wo) := r{z, W{z, Wo)). 

As compared to the standard ODE construction of, e.g., [B, VI, TZl, Z5], 
in which (2.8)-(2.9) together arc considered as a fixed-point equation for the 
joint variable {w, z), this amounts to treating w implicitly. This is a standard 
device in situations of limited regularity; see, e.g., [CP, GMWZ, RZ]. 

It remains to show, first, that W, hence T, is well-defined on a space of 
slowly-exponentially-growing functions and, second, that T is contractive on 
that space, determining a C'^ solution z = z{wq) similarly as in the usual 
ODE construction. We carry out these steps in the following subsections. 
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2.3 Nonlinear energy estimates 

Define now the negatively-weighted sup norm 

||/||_, :=supe-^*|/(t)|^M, 
t>o 

noting that \f (1)1^2 < e^*||/||_0 for all t > 0, and denote by B^^f the Banach 
space of functions bounded in || • norm. Define also the auxiliary norm 

ll/IU. :=supe-''*|/(i)|^.. 

Lemma 2.4 ([Z5]). Under assumptions (Al)-(A3), (H0)-(H2), for all 1 < 
p < oo, < r < 4, 

(2.12) |n„|ip_>vF'-'P , \^cs\wr,p-*W'^'P ^ C ■ 

Proof. Recalling that L has at most finitely many unstable eigenvalues, we 
find that may be expressed as 

p 

where (f)j, j = are generalized right eigenfunctions of L associated 

with unstable eigenvalues Aj, satisfying the generalized eigenvalue equation 
{L — \jY^(f)j = 0, Tj > 1, and (f)j arc generalized left eigenfunctions. Noting 
that 0j, (pj and derivatives decay exponentially by standard theory [He, ZH, 
MaZl], and estimating 

\diUJ\LP = \J2^JAy)\LP < E \dicl>j\LSMf\LP < C\f\LP 
j j 

for 1/p + 1/q = 1, we obtain the claimed bounds on n„, from which the 
bounds on lies = Id — n„ follow immediately. □ 

Proposition 2.5. Assuming (Al)-(A3), (H0)-(H2), denote A^ := A^{u{x)) 
and take without loss of generality a{Aii < 0. Then, there exist a skew- 
symmetric matrix K{x) and constants C > 0, C^, > such that, defining the 
scalar weight a by a{0) = 1 and = C^,\Ux\a, the quadratic form 

S{w) ((C^^o + CKd^ - A'^dDaw, w) 



11 



is equivalent to and, for Wt — Lw 



(2.13) dt£{w) < -eiwlj,,,, + C(|/|i2 + \w\l,). 

Proof. This follows by a linear version of the argument for Proposition 4.15, 
[Z3], somewhat simplified by the stronger block structure assumptions (Al)- 
(A3) made here. We first observe that, by symmetry of {A'^A)ii and {A'^A)±, 

A'A(x)-A(x) + ( ° ^(l^^l)^ 
A A{x) - A{x) + ^^^^^^^^ 0(|[/.|)J 

where A is the symmetric part of A^A. Next, we recall from [Kaw, KSh] that 
(A2)-(A3) imply existence of skew-symmetric K± such that ^{KA+B)± > 0. 
Defining K{x) as a smooth interpolant between Kj- defined by arc length 
along the path of U between U±, we obtain the result by a straighforward (if 
somewhat lengthy) computation for C > and C* >> C C > sufficiently 
large, where "good" terms 

-{dlw,aA°Bdlw >< -^(a>2, aa>2|, 

-C^{d^w,aA^Bd^w >< -C^9{d^W2,ad^W2\, 
-C{d^w,^{KA + A^B)ad^w) < -Ce{d^w,ad^w) 

+ C2{d^wi, \tJ^\ad^wi) + C2{d^W2, ad^W2), 

and 

-{d^wi, {aA^A)-^idlwi) = {d^wi,d,^{aA^A)ii ■ d^wi) 

< -6C^{dxWi, \ Uj:\adxWi) 

resulting by rearrangement /integration by parts of the various components 
of dtS{w), when summed together, after absorbing remaining terms, give a 
contribution of —6*1^1^1,2 + Cd^l^g + 1/1^2), verifying the result. See [MaZ4, 
Z2, Z3] for further details. □ 

Corollary 2.6. Under assumptions (A1)-(A3), (H0)-(H2),for\\z\\i2 bounded 
and 6 sufficiently small, the solution w of (2.10) exists for allt>0 and, for 
any constant 9 > and some C = C{6), 
(2.14) 



e-''\w\l,,2{s) ds + e-''\wt\jii < C\wo\l. + C I e-''{\w\l, + \z\l2){s) ds- 

Jo 
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likewise, for 9 > 2r], 
(2.15) 

/+00 r+oo 
e'^'-'^\w\%r4s) ds < C\w\%,{t) +C e'^'-'\\w\l, + \z\l,){s) ds 

<Ce''\wo\li+Ce^'^\\\w\\l, + \\z\\l, ). 



Proof. Observing that 







(2.16) \U,,N\w + z)\l2 < 

we obtain by Proposition 2.5 the inequaUty 

dtS{w) < -9\w\hi,2 + CS{\w\jji,2 + \z\jji,2) + C{\w\l2 + \z\l2) 
for some 9 > 0, which, for 6 sufficiently small, gives 

dtS{w) < -^1^01^1.2 + C2{\w\l2 + 1^112) 

by |-2|ifi,2| < C\z\l2 (equivalence of finite-dimensional norms). From this in- 
equality, global existence follows immediately and (2.14) and the first hue 
of (2.15) follow readily by Gronwall's inequality together with equivalence of 
S and lit'll,!. For further discussion, see [Z6], proof of Proposition 2.6. 

Bounding |w(t)|l^i in the second line of (2.15) using (2.14) now yields 
(2.17) 

/+OO f+OO 
e'^'-'^\w\l,,2{s)ds<Ce''(\wo\]ii + J e-''{\w\l2 + \z\l2){s) ds) , 

whereupon the final hue of (2.15) then follows by 

r+oo , r+oo 

/ e'^'-'\\w\l2 + \z\l2){s)ds<{e'' / e^^""^)^ ds) (||w||2.^ + J 

<Ce'^\\\w\\,._^ + \\z\\,.J. 

□ 

Corollary 2.7. Under (A1)-(A3), (H0)-(H2), for \\zi\\li^, W^^hl^ bounded 
and 5 sufficiently small, solutions wi, zi andw2,Z2 of (2.10) exist for allt > 
and, for any constant 9 > and some C — C{9), 



(2.18) 



/ e I wi - W2 1 1,1,2(5) (is < C I Wo,i - Wo,2 1 Ifi 
^0 

+ C e~^'{\wi - W2\l2 + \zi - Z2\l2){s) ds; 
Jo 
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likewise, for 9 > 2r), 

/+00 
e^(t-')\wi - W2\li,2{s) ds < C\wi - W2\li{t) 

/ + 00 
e^^^-'\\wi - W2\l2 + \zi - Z2\l2){s) ds 

< Ce^* 1^0,1 - wo,2\m 
+ Ce''''\\\wi - W2\\l2 +\\zi-Z2\\l2 ). 

-V -V 

Proof. Subtracting the equations for wi^^zi and W2,2;2, we obtain, denoting 
w :— wi — W2, z :— zi — Z2, the equation 

(2.20) Wt-Lw = U,,{N\wi + zi) - N\w2 + Z2)) 

with initial data wq = Wo,i ~ ^0,2 at t = 0, whence the resuh follows by the 
Lipshitz bound 6 on A^'' : H^'^ — > and the same argument used in the 
proof of Corollary 2.6. □ 

Corollary 2.8. Under assumptions (A1)-(A3), (H0)~(H2), for3uj<r)<(3 
and 5 > and Wq G sufficiently small, for each z e B^r), there exists a 
unique solution w —: 'W{z,Wo) e of (2.9), (2.10), with 

(2.21) \MlI^ < C{\wo\m,2 + 6\\z\\_rt) 
and 

(2.22) \\Wizi,Wo,l)-Wiz2,Wo,2)\\Ll^ < C5\\zi-Z2\\-r, + C\\wo,l-Wo,2\\m,2. 

Proof. We have already shown existence, while uniqueness follows from (2.22). 
Thus, we need only verify the estimates. 

Consider a pair of data zi,wo^i and ,22,^0,2, and compare the resulting 

solutions, denoting {z,w,Wq) := {zi — Z2,Wi — W2,Wq^i — Wq^2)- Using the 
integral representation (2.9), and applying (2.3)(i), (2.6)(i), and the definition 
of II • II we obtain for alH > 

(2.23) \w(t)\L2 < Ce'^l*l|wo|L2 + C5 [ e'"l*-^l(|^i;jj^2 + \z\h2){s) ds. 

Jo 

Estimating 

C6 f e'^^'-'^\z\H2{s)ds <C6\\z\\_rj f e^^'-^^e"^ cis < C5||i||_^e''* 
Jo Jo 
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and, by (2.18) with 9 — Sou together with the Cauchy-Schwarz inequahty, 

CS £ e-(*-^) \w\„2 is) ds<CS(^ £ e-'^(*-^)) ( e'^^'-'^ \w\l. {s) ds) ^'^ 

we obtain, substituting in (2.23), 
(2.24) 

\w{t)\l. < Ce'-^'\\wo\l. + C6'[\\4\e'^' + f e'-'^'-'WAl. + \z\l.){s)ds) 

<e^^*(cKl^i + C5^(||^|li.^ + ||i|r-.)). 
This yields 

\\w{t)\\L2_^ < C\wo\^^ + C5{\\w\\li^ + ||i||_^), 

from which (2.22) follows by smallness of 5. The bound (2.21) follows simi- 
larly. □ 

2.4 Basic existence result 

Proof of Proposition 1.2 without translational independence, (i) (Contraction 
mapping argument) We find using (2.22), (2.6)(ii), and Lemma 2.2 that 

11^(^1, ^i^O,l) - f{z2,Wo,2)\\-n 

/+00 
e'^^'-'\\w\H^ + \z\h^){s) ds 

< supCi(5(^||i||_^ + e-''* r e^(*-*)|w|i/2(s)fis). 
Using the Cauchy-Schwarz inequality and (2.19) with 9 = P to estimate 

"+00 . p+oo ^ 1 /9 ^ r+oo 



< Cs{\w\Ut) + j\^^'-^\\w\h + \z\l.){s)ds)''" 

< Cae'^M |wo|//i + ||w||l2 + ||i;||L2 ) 
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and applying (2.22), we obtain 

(2.26) ||t(2;i,tt;o,i) - t(z2, wo,2)||-„ < C6{\\wo\\m,2 + \\z\\lI^)- 
Parallel estimates yield 

\\f{z,wom\\-r, < CS{\wo\m,2 + S\\z\\^,), 

so that, taking 5 and jwolffi.a sufficiently small, that T{-,wo) maps the ball 
B{0,r) C B-rj to itself, for r > arbitrarily small but fixed. 

This yields at once contractivity on 5(0, r), hence existence of a unique 
fixed point z = Z{wq), and Lipshitz continuity of Z from to by the 
Banach Fixed- Point Theorem with Lipschitz dependence on parameter wq. 

(a) (Existence of a Lipschitz invariant manifold) Defining 

r+oo 

(2.27) $(wo) := 2(wo)|t=o = - / e-^^n„iV^(^;(s)) ds, 

Jo 

we obtain a Lipschitz function from E„, whose graph over 5(0, r) 

is the invariant manifold of solutions of (2.7) growing at exponential rate 
\v{t)\ < Ce^* in forward time. Prom the latter characterization, we obtain 
evidently invariance in forward time. Since the truncated equations (2.7) 
agree with the original PDE so long as solutions remain small in H^'^, this 
gives local invariance with respect to (1.2) as well. By uniqueness of fixed 
point solutions, we have Z(0) — and thus $cs(0) = 0, so that the invariant 
manifold passes through the origin. Likewise, any bounded, sufficiently small 
solution of (2.1) in H^'^ is a bounded small solution of (2.7) as well, so by 
uniqueness is contained in the center stable manifold. 

(zzz) (Quadratic-order tangency) By (2.27), (2.3), (2.4), (2.15), amd (2.21), 

r+oo 

Mwcs)\m,2^ / e-'^'U^N^w + z)(s) ds 

r+oo 

<C / e-^\\W\]ju2 + \Z\]j,,,){s)ds 
Jo 

<C2{\Wcs\li,2 + \\Zf_^). 

By Z{0) — and Lipshitz continuity of Z, we have < C\wcs\m-^, 

whence 

(2.28) |$(wes)|i^M < CsIu^csIhm, 

verifying quadratic-order tangency at the origin. □ 
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2.5 Translation-invariance 

We conclude by indicating briefly how to recover translation-invariance of 
the center stable manifold, following [Z5, TZl]. Differentiating with respect 

to X the traveling- wave ODE, we recover the standard fact that := t/j; is 
an zero eigenfunction of L, by the decay of Ux noted in Remark 1.5. 
Define orthogonal projections 

(2.29) n2:=4^, ni:=id-n2, 

onto the range of right zero-eigenf unction (f) :— Ux oi L and its orthogonal 
complement 0-*- in L^, where (•, •) denotes standard inner product. 

Lemma 2.9. Under the assumed regularity h e C'''^^, k > 2, Uj, j — 1,2 
are bounded as operators from to itself for Q < s <k -\-2. 

Proof. Immediate, by the assumed decay oi ()) = Ux and derivatives. □ 

Proof of Proposition 1.2, with translational independence. Introducing the shifted 
perturbation variable 

(2.30) V{x,t) ■.^U{x + a{t),t) -U{x) 
we obtain the modified nonlinear perturbation equation 

(2.31) dtV ^LV +N{y)x-dta{(t> + dxV), 

where L :— ^{U) and M as in (2.2) is a quadratic-order Taylor remainder. 

Choosing dta so as to cancel 112 of the righthand side of (2.31), we obtain 
the reduced equations 

(2.32) dtV = n^{LV + N{V)) 
and 

t:2{LV + N{V)) 



(2.33) dta = 



l + T:^{dxV) 



V e 0-*-, 7r2F := (0, y) 1 01 ^2, of the same regularity as the original equations. 
y(o) e 0^, or 

(0,C/o(x + a)-t/(a;)) = 0. 
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Assuming that Uq lies in a sufficiently small tube about the set of translates 
of U, or Uo{x) = (U + W){x — j3) with |iy|j:/i,2 sufficiently small, this can be 
done in a unique way such that a :— a — P is small, as determined implicitly 



With this choice, translation invariance under our construction is clear, with 
translation corresponding to a constant shift in a that is preserved for all 
time. 

Clearly, (2.33) is well-defined so long as |9xV^|l2 < C|y|i^i,2 remains 
small, hence we may solve the v equation independently of a, determining 
a-behavior afterward to determine the full solution 



Moreover, it is easily seen (by the block-triangular structure of L with re- 
spect to this decomposition) that the linear part IIiL = riiLrii of the in- 
equation possesses all spectrum of L apart from the zero eigenvalue asso- 
ciated with eigenfunction 0. Thus, we have effectively projected out this 
zero-eigenfunction, and with it the group symmetry of translation. 

We may therefore construct the center stable manifold for the reduced 
equation (2.32), automatically obtaining translation-invariance when we ex- 
tend to the full evolution using (2.33). See [TZl, Z5] for further details. □ 

3 Conditional stability analysis 

Define similarly as in Section 2.5 the perturbation variable 



for U a solution of (1.2), where a is to be specified later. Subtracting the 
equations for U{x + a{t),t) and U{x), we obtain the nonlinear perturbation 
equation 




9a,e?(O,O) = (0,a,C7(-)) = 10112^0. 



U{x, t) = U{x - a{t)) + V{x - a{t),t). 



(3.1) 



V{x, t) := U{x + a{t),t) - U{x) 



(3.2) 



Vt-LV =Af{V),-a{U, + V,), 
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where L := —Ad^ + d^Bd^ as in (1.6) denotes the linearized operator about 
U and jV as in (2.2) is a nonlinear residual, satisfying 

.33. Af{V)^0{\V\' + \Vm). 

^ ■ ^ d,J\f{V)^0{\V\\d,V\ + \V,\'' + \V\\V,,\) 

so long as (hence \V\l'=° and |t/|L°°) remains bounded. 



3.1 Projector bounds 

Let Uu denote the eigenprojection of L onto its unstable subspace S„, and 
lies = Id — n„ the eigenprojection onto its center stable subspace T,cs- 

Lemma 3.1. Assuming (HO)-(Hl), 
(3.4) n,9, = 

for j = u, cs and, for all 1 < p < 00, < r < 4, 

|n„|vF'-,p_>vK'-,p, |nu|vK'-,p_»v^r-,p < c, 

(3-5) ~ ~ 

Proof. Recalling (see the proof of Theorem 1.2) that L has at most finitely 
many unstable eigenvalues, we find that 11^ may be expressed as 

where (f)j, j = l,...p are gcncrahzed right eigenfunctions of L associated 
with unstable eigenvalues satisfying the generalized eigenvalue equation 
(L — XjY^4>j = 0, rj > 1, and (pj are generalized left eigenfunctions. Noting 
that L is divergence form, and that Xj 7^ 0, we may integrate {L — XjY^cjyj — 
over R to obtain A^-^ / (f)jdx — and thus / (f)jdx — 0. Noting that 0^, 4>j and 
derivatives decay exponentially by standard theory [He, ZH, MaZl], we find 
that 

(j). = d^^j 

with $j and derivatives exponentially decaying, hence 

3 



19 



Estimating 



j j 

for l/p+ 1/q = 1 and similarly for d^Iluf, we obtain the claimed bounds on 
and n„, from which the bounds on Ucs = Id — n„ and lies = W — n„ 
follow immediately. □ 

3.2 Linear estimates 

Let Gcs{x, t; y) := Ilcse^^Sy{x) denote the Green kernel of the linearized solu- 
tion operator on the center stable subspace J^cs- Then, we have the following 
detailed pointwise bounds established in [TZ2, MaZl]. 

Proposition 3.2 ([TZ2, MaZl]). Under (A0)-(A3), (H0)-(H3), (D1)-(D3), 
the center stable Green function may he decomposed as Gcs — E + H + G, 
where 

(3.6) E{x,t;y) = d,Uix)e,iy,t), 

(3.7) eiy,t) = ^ ^errfn (^^) - errfn (^^)) l^iv) 

for y < and symmetrically for y > 0, e R" constant, 

J 

(3.8) H(x,t;y) = J20{e-''>')6,--ap{-y), 9 > 0, 
and 

n 

\G{x,t;y)\ < Ce-^^l^-^/l+t) +^^-l/2g-(x-j/-a^t)VMtg-r,a.+ 

fe=l 

(3.9) ' ^ 



I \^ V .-1/2 -ix-a-it-\y/a^ my Mt -nx^ 



>0, aj <0 



I \^ V f-1/2 -{x-a 



(t-\v/H\)f/Mt -nx- 

'^i\%t\^\y\r 

a^>0,at>0 
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(3.10) 

n 
k=l 

I V V f-^/2p-{x-a-{t-\y/a^\))yMt -VX+ 

+ 2^ ^{Kt\>\y\r e ^ e 



aj, >0, a^. <0 



I V V ^-1/2 -(x-at(t-|y/a^|))VMt -„a;-\ 



a^. >0,at>0 



/or y < and symmetrically for y > 0, for some rj, C , M > 0, where are 
the eigenvalues of Aj^ = dF{U±), a*,j = 1, . . . , J are the eigenvalues of An, 
denotes the positive/negative part of x, and indicator function X{\a~t\>\y\} 
is 1 for |a^t| > \y\ and otherwise. 



Proof. As observed in [TZ2], it is equivalent to establish decomposition 

(3.11) G^G^ + E + H + G 

for the full Green function G{x,t;y) := e^^6y{x), where 

p 

Gu{x,t;y) := n„e%(x) = e^* J] </.,(a;)0,(|/)* 

for some constant matrix M G C^^^ denotes the Green kernel of the linearized 
solution operator on E„, (pj and (pj right and left generalized eigenfunctions 
associated with unstable eigenvalues Xj, j — 1, . . . ,p. 

The problem of describing the full Green function has been treated in 
[ZH, MaZ3] , starting with the Inverse Laplace Transform representation 

(3.12) G{x, t; y) = e^^dy{x) = j e^\\ - L{e))-^5y{x)d\ , 
where 

r := d{X : 5RA < r^i - 772|$^A|} 

is an appropriate sectorial contour, r/i, 772 > 0; estimating the resolvent kernel 
G\{x, y) :— (A — L{e))''^5y{x) using Taylor expansion in A, asymptotic ODE 
techniques in ,x, y, and judicious decomposition into various scattering, ex- 
cited, and residual modes; then, finally, estimating the contribution of various 
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modes to (3.12) by Riemann saddlepoint (Stationary Phase) method, moving 
contour F to a optimal, "minimax" positions for each mode, depending on 
the values of {x,y,t). 

In the present case, we may first move F to a contour F' enclosing (to the 
left) all spectra of L except for the p unstable eigenvalues Xj, j — 1, . . . ,p, 
to obtain 

G{x, t;y)= I e^*(A - L)-^d\ + J] ResidueA^.(e) (e^*(A - L)-^5y{x)) , 

where ResidueA^(£) (e'^*(A — L)^^8y{x)) — Gu{x,t;y), then estimate the re- 
mammg term fp, e^\X- L)-^dX on minimax contours as just described. See 
the proof of Proposition 7.1, [MaZ3], for a detailed discussion of minimax es- 
timates E+G and of Proposition 7.7, [MaZ3] for a complementary discussion 
of residues incurred at eigenvalues in {3?A > 0} \ {0}. See also [TZl]. □ 

Corollary 3.3 ([MaZl]). Assummg (Al)^(A3), (H0)-(H3), (D1)-(D3), 

/ + 00 
H{;t-y)f{y)dy\L.<Ce-'y\L., 
-oo 



/+00 ^ ^ ^ 

G{-,t-y)f{y)dy\L.<C{l+t)--^-^-^\f\L., 
-oo 



(3.15) I / G,(-, ^; y)f{y)dyy < C{1 + i)-i(i4)-i 

J —oo 

for all t > 0, some C > 0, for any 1 < q < p ( equivalently, 1 < r < p) and 
f e LI, where 1/r + 1/q ^ 1 + 1/p. 

Proof. Standard convolution inequalities together with bounds (3.8)-(3.10); 
see [MaZl, MaZ2, MaZ3, Z2] for further details. □ 

Corollary 3.4 ([Z4]). The kernel e satisfies 

\e,i;t)\L.,\eti;t)\Lr<Cr^^'-'/^\ 
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for all t > 0. Moreover, for y <0 we have the pointwise bounds 

_1 ^-^ / (i; + a~f)^ fa-a~f)2 X 

\ey{y,t)\,\et{y,t)\ <Ct ^ 2^ (e Mt +e m* j, 



for M > sufficiently large, and symmetrically for y > 0. 

Proof. Direct computation using definition (3.7); see [Z4, MaZ2, MaZ3] or 
[Z5], Appendix A. □ 

3.3 Reduced equations II 

RecaUing that dj U is a stationary solution of the hnearized equations Ut = 
LU, so that LdxU = 0, or 

/oo 
G{x,t;y)UMdy^e'^'U,{x)^dM{x), 
-oo 

we have, applying Duhamel's principle to (3.2), 

/oo 
G{x,t;y)Vo{y) dy 
-OO 

/•t poo 

- Gy{x,t - s]y){M{V) + aV){y,s)dyds + a{t)dxu{x). 

Jo J-oo 

Defining 

/oo 
(^{v,'t)yo{.y) dy 

+ ey{y,t-s){M{V) + aV){y,s)dyds, 

Jo J-oo 

following [ZH, Z4, MaZ2, MaZ3], where e is as in (3.7), and recalling the 
decomposition G — E+H+Gu+G of (3.11), we obtain the reduced equations 

/oo 
{Gu + H + G){x,t;y)Vo{y) dy 
-OO 

nt POO 

(3.17) -/ / {Gu + G)y{x,t-s;y){Af{V) + aV){y,s)dyds 

Jo J-oo 
rt roo 

+ 11 H{x,t- s;y){M{V)+aV)y{y,s)dyds, 

Jo J-oo 
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and, difTerentiating (3.16) with respect to t, and observing that ey{y,s) —r 
as s — > 0, as the difference of approaching heat kernels, 

/oo 

+ / / eyt{y^- s){M{y) + aV){y,s)dyds. 

Jo J-oo 



Note that this (nonlocal in time) choice of a and the resulting reduced 
equations are different from those of Section 2.5. As discussed further in 
[Go, Z4, MaZ2, MaZ3, Z2], a may be considered in the present context as 
defining a notion of approximate shock location. 



3.4 Nonlinear damping estimate 

Proposition 3.5 ([MaZ3]). Assuming (A0)-(A3), (H0)-(H2), let Vq E H^' , 
and suppose that for < t < T, the norm of V remains hounded by a 
sufficiently small constant, forV as in (3.1) andu a solution of (1.2). Then, 
for some constants 6*1,2 > 0, for all < t < T , 

(3.19) \\Vml^<Ce-'^'\\Vml^ + C fe-'^^'-'\\V\l. + \a\'){s)ds. 

Jo 

Proof. Essentially identical to the proof of Proposition 2.5 and corollaries, 
but using Moser's inequality to bound the nonlinear term instead of the 
Lipshitz bound imposed by truncation. See [MaZ4, Z2, Z3] for detailed proofs 
of more general results. □ 

3.5 Proof of nonlinear stability 

Decompose now the nonlinear perturbation V as 

(3.20) V{x,t) =w{x,t) + z{x,t), 
where 

(3.21) w := U,sV, z := U^V. 
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Applying Ucs to (3.17) and recalling commutator relation (3.4), we obtain 
an equation 



w{x,t) = G{x,t;y)wo{y)dy 

J —OO 

rt roo 

(3.22) - / / Gy{x,t - s;y)IlcsiXiV) + aV){y,s)dyds 

Jo J -oo 

H(x, t - s; y)Ucs(J^(V) + aV)y(y, s)dy ds 



oo 

t POO 



10 J-oo 

for the flow along the center stable manifold, parametrized hy w & T^cs- 

Lemma 3.6. Assuming (Al)-(A3), (H0)-(H2), forV lying initially on the 
center stable manifold Ai cs, 

(3.23) l-^lw'P < C'|u'|^2 

for some C > 0, for all 1 < p < oo and < r < 4, so long as \w\h^ remains 
sufficiently small. 

Proof By (1.9), we have immediately |2;|_ffi,2 < 6*1^1^1,2, whence (3.23) fol- 
lows by equivalence of norms for finite-dimensional vector spaces, applied 
to the p-dimensional subspace S„. (Alternatively, we may see this by di- 
rect computation using the explicit description of HuV afforded by Lemma 
3.1.) □ 

Proof of Theorem 1.3. Recalling by Theorem 1.2 that solutions remaining for 
all time in a sufficiently small radius neighborhood M of the set of translates 
of u lie in the center stable manifold A^cs, we obtain trivially that solutions 
not originating in Ales must exit H in finite time, verifying the final assertion 
of orbital instability with respect to perturbations not in M.cs- 

Consider now a solution V G A^cs, or, equivalently, a solution w e Ecs of 
(3.22) with z = ^cs{w) G S„. Define 

(3.24) C(^) — sup f|w|^j3(l + s)i + (|«;|ioc + |Q;(s)|)(l + s)5). 

We shall establish: 

Claim. For all t > for which a solution exists with C, uniformly bounded 
by some fixed, sufficiently small constant, there holds 

(3.25) C{t)<C2{E^ + C{tf) for -.^ \VoW^h^ . 
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Prom this result, provided Eq < 1/4C|, we have that ({t) < 2C2E0 
imphes ({t) < 2C2E0, and so we may conclude by continuous induction that 



(3.26) 



C{t) < 2C2E0 



for all t > 0, from which we readily obtain the stated bounds. (By standard 
short-time H'' existence theory, V G exists and ( remains continuous so 
long as C remains bounded by some uniform constant, hence (3.26) is an open 
condition.) 

Proof of Claim. By Lemma 3.6, 



whence 



\wo\L^nHS < CEq 



Likewise, by Lemma 3.6, (3.24), (3.3), and Lemma 3.1, for < s < i and 
2 <p < 00, 



(3.27) 



\UUNiV) + aV)iy, s)\l2 < CC(t)'(l + s)^^ 
\U,,{N{V) + aV)y{y, s)\l. < CC(t)'(l + s)"^ , 



Combining the latter bounds with representations (3.22)-(3.18) and ap- 
plying Corollary 3.3, we obtain 
(3.28) 



\w{x,t)\LP < 



G{x,t;y)wQ{y) dy 



LP 



t POD 



+ 
+ 

< Eo{l + t) 



J -00 

t POD 



J -00 



Gy{x, t - s; y)ncs{N{V) + aV){y, s)dy ds 
H{x, t - s; y)Ucs{N{V) + aV)y{y, s)dy ds 

rt 

+ CC{tf / {t- s)~~^^^{l + s)--^dyds 
Jo 



LP 



LP 



-e{t-s) , 



<C{Eo + aty){l + t)-^''-p 



1 + s) 2 dy ds 
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and, similarly, using Holder's inequality and applying Corollary 3.4, 
(3.29) 

/oo 
\et{y,t)\\Vo{y)\dy 
-oo 

pt P+OO 

+ \eyt{y,t-s)\\fl,,{N{V) + aV){y,s)\dyds 

Jo J -oo 

< letlL-miLi + CCitf f \eyt\L2{t - s)\Il,,{N{V) + aV)\L2{s)ds 

Jo 

<Eoil + t)--^+CC{tf [ + 

Applying Lemma 3.5 and using (3.28) and (3.29), we obtain, finally, 
(3.30) \w\Hs{t) < C{Eo + C{tf){l + t)-l 

Combining (3.28), (3.29), and (3.30), we obtain (3.25) as claimed. 

As discussed earlier, from (3.25), we obtain by continuous induction 
(3.26), or ^ < 2C2I Vo|Lin//2, whereupon the claimed bounds on \V\lp and 
\V\h3 follow by (3.28) and (3.30), and on \a\ by (3.29). Finally, a computa- 
tion parallel to (3.29) (see, e.g., [MaZ3, Z2]) yields \a{t)\ < C{Eo + C{ty), 
from which we obtain the last remaining bound on |Q;(i)|. □ 

Acknowledgement. Thanks to Milena Stanislavova and Charles Li for 
two interesting discussions that inspired this work. 

A Example systems 

1. The general Navier-Stokes equations of compressible gas dynamics, writ- 
ten in Lagrangian coordinates, appear as 

{Vt-Ua: = 0, 
(e + u^/2)t + (pu)^ = {{k/v)T^ + {ii/v)uu^)^, 

where v > denotes specific volume, u velocity, e > internal energy, 
T — T{v, e) > temperature, p — p{v, e) pressure, and /i > and k > are 
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coefficients of viscosity and heat conduction, respectively. Defining f/i := (v), 
U2 := (M,e+ |?i|V2), we find that conditions (A1)-(A3) and (H1)-(H3) are 
thus satisfied under the mild assumptions of ideal temperature dependence 

(A.2) T = T{e) 

monotone temperature-dependence 

(A.3) Te > 0, 

and thermodynamic stability of the endstates, 

(A.4) (p,)± < 0, (7;)± > 0; 

see, e.g., [MaZ4, Z2, Z3, TZ3] for further discussion. Notably, this allows the 
interesting case of a van der Waals-type equation of state, with p„ > for 
some values of v along the connecting profile. 

2. The equations of MHD in Lagrangian coordinates are 
(A.5) 

vt - Mix = 0, 

uit + (P + (l/2/xo)(5| + Bl))^ = {{u/v)ui^)^, 

U2t - {{1/ lXo)BlB2)x = {{l^/v)u2:,)x, 
U3t - {{l/^lo)B*B3)x = ((z//f )m3x)x-, 

{vB2)t - {Blu2)x = iil/(rf^ov)B2x)x, 

(e + (1/2)K + ul + ul) + {l/2fio)v{Bi + Bl)\ 

+ + (l/2/io)(5| + Bl))u^ - {l/pio)B*{B2U2 + BsU3)l 

= [{u/v)uiUix + {l^/v){u2U2x + U^U^x) 

+ {k/v)Tx + {l/aiJ,lv){B2B2x + B3Bsx)]x, 

where v denotes specific volume, u — (1^1,1*2,1*3) velocity, p — P{v,e) pres- 
sure, B = {B^, B2, B3) magnetic induction, constant, e internal energy, 
T = T{v, e) > temperature, and /x > and u > the two coefficients 
of viscosity, k. > the coefficient of heat conduction, /xq > the magnetic 
permeabihty, and a > the electrical resistivity. Under (A.2)-(A.4), condi- 
tions (A1)-(A3) are again satisfied, and conditions (H1)-(H3) are satisfied 
(sec [MaZ4]) under the generically satisfied assumptions that the endstates 
be strictly hyperbolic (i.e., have simple eigenvalues), and the speed s be 
nonzero, i.e., the shock move with nonzero speed relative to the background 
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fluid velocity, with Ui := (v), U2 := {u,B,e+ \u\'^/2 + v\B\'^/2iio)- (For gas 
dynamics, only Lax-type shocks and nonzero speeds can occur, and all points 
U are strictly hyperbolic.) 

3. (MHD with infinite resistivity/permeability) An interesting variation 
of (A. 5) that is of interest in certain astrophysical parameter regimes is the 
limit in which either electrical resistivity a, magnetic permeability /Xq, or 
both, go to infinity, in which case the righthand sides of the fifth and sixth 
equations of (A. 5) go to zero and there is a three-dimensional set of hyperbolic 
modes {v, VB2, vB^) instead of the usual one. By inspection, the associated 
equations are still linear in the conservative variables. Likewise, (Al)-(A3), 
(H1)-(H3) hold under (A.2)-(A.4) for nonzero speed shocks with strictly 
hyperbolic endstates. 

4. (multi-species gas dynamics or MHD) Another simple example for 
which the hyperbolic modes are vectorial is the case of miscible, multi-species 
flow, neglecting species diffusion, in either gas dynamics or magnetohydro- 
dynamics. In this case, the hyperbolic modes consist of k copies of the 
hyperbolic modes for a single species, where k is the number of total species. 
Again, (Al)-(A3), (H1)-(H3) hold for nonzero speed shocks with strictly 
hyperbolic endstates under assumptions (A.2)-(A.4). 
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